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Abstract 

We developed a three-dimensional gaseous tracking device and performed a direction- 
sensitive dark matter search in a surface laboratory. By using 150 Torr carbon- 
tetrafluoride (CF4) gas, we obtained a sky map drawn with the recoil directions of 
the carbon and fluorine nuclei, and set the first limit on the spin-dependent WIMP 
(Weakly Interacting Massive Particles)-proton cross section by a direction-sensitive 
method. Thus, we showed that a WIMP-search experiment with a gaseous track- 
ing device can actually set limits. Furthermore, we demonstrated that this method 
will potentially play a certain role in revealing the nature of dark matter when a 
low-background large-volume detector is developed. 

Key words: Time projection chamber, Micro pattern detector, dark matter, 

WIMP, direction-sensitive 

PAGS: , 14.80.Ly 29.40.Cs, 29.40.Gx, 95.35. +d 



Email address: miuchi@cr.scphys.kyoto-u.ac.jp (Kentaro Miuchi"). 
Preprint submitted to Elsevier Science 1 February 2008 



1 Introduction 



Revealing the nature of dark matter has received much more attention since 
the successive results of the WMAP cosmic microwave background all-sky 
observation [1]. Many dark matter direct search experiments have been concluded[2,3,4,5,6,7], 
are currently being performed[8,9,10,ll,12,13,14], or are in the planning stage[15]. 
Weakly interacting massive particles (WIMPs) are very plausible candidates 
of dark matter. Most dark matter search experiments are designed to measure 
only the energy transferred to the nucleus through WIMP-nucleus scatterings. 
In these experiments, the distinct signals of WIMPs are the annual modula- 
tion of the energy spectrum, the halo is isotropic. Because the amplitude of 
an annual modulation signal is only a few % of the rate, positive signatures 
of WIMPs are very difficult to detect with these detectors [4]. On the other 
hand, the motion of the solar system with respect to the galactic halo is con- 
sidered to provide a much larger asymmetry in the direction-distribution of 
the WIMP velocity observed at Earth[16]. A detection method of a positive 
signature of WIMPs by measuring the directions of the nuclear recoils have 
been studied[2,6,14,17]. If the distribution of the local halo were not isotropic 
- perhaps quite clumpy- direction-sensitive methods would then potentially 
provide information on it. 

Gaseous detectors are one of the most appropriate devices to detect this 
"WIMP-wind" [18,19]. The DRIFT project has performed underground mea- 
surements for more than two years with 1 m^ time projection chambers (TPC) 
filled with a low-pressure CS2 gas aiming to detect the mainly WIMP-wind via 
spin-independent (SI) interactions [14]. Since former studies suggest to perform 
the WIMP-search both via SI and spin-dependent (SD) interactions [20, 21], we 
proposed a SD-sensitive WIMP-wind search, named NEWAGE (NEw genera- 
tion WIMP-search with an Advanced Gaseous tracking device Experiment) [22] . 
We developed a /i-TPC with a detection volume of 23 x 28 x 30 cm^ and 
performed a direction-sensitive WIMP-search experiment in a surface labora- 
tory. 



2 Detector performance 

In this section, we describe the detector performance, which is also described in 
Ref . [23] . We mainly measured the detector response in the energy range that 

we used for the WIMP-search analysis (100-400keV, DM energy range). The 
lower threshold (lOOkeV) was set at the lowest energy for an effective gamma- 
ray rejection with a 150 Torr CF4 gas, and the higher threshold (400keV) was 
set at the expected highest recoil energy by the WIMPs at the escape velocity. 
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2.1 ii-TPC 



A //-TPC is a gaseous time projection chamber with a micro pixel chamber (//- 
PIC) readout. A /x-PIC is a gaseous two-dimensional position-sensitive detec- 
tor manufactured by printed circuit board (PCB) technology. PCB technology 
realizes an economical mass production, which is one of the most important 
requirements to fabricate a dark matter detector. After a general-purposed de- 
tector R&D, we studied the performance of a /x-TPC with a detection volume 
of 10 X 10 X lOcm^ using a 150 Torr CF4 gas as our first step towards a dark 
matter experiment [24]. In this study, we demonstrated the tracking of nuclear 
recoils and the gamma-ray rejection with a CF4 gas of 150 Torr. We also found 
that we need some more studies to operate the /x-TPC with a CF4 gas at a de- 
sign value of 30 Torr for more than a few days. We then developed a large-sized 
(31 X 31 cm^) /x-PICs[25], and subsequently developed a large- volume //-TPC 
with a detection volume of 23 x 28 x 30 cm"'. We studied the fundamental 
properties of this large-volume /x-TPC with a standard Ar-C2H6 gas mixture 
at normal pressure[26]. The size is of the /x-PIC is restricted by the working 
size of the PCB, and we plan to realize a larger readout area by placing many 
pieces, for instance 4x4 pieces for a 1 m^, of /i-PICs. Before we take the 
next technology steps, we performed a first direction-sensitive spin-dependent 
dark matter search experiment to demonstrate that a gaseous detector and a 
direction-sensitive method can actually set limits on the WlMPs. We report 
on the results of a direction-sensitive dark matter search experiment using 
this large-volume //-TPC with a 150 Torr CF4 gas (9.0g effective mass) in a 
surface laboratory. 

We used the /x-TPC described in Ref. [26], and thus most parts of the detector 
were unchanged. Herein, we describe only the significant properties and the 
modified parts of the detector system. Figure 1 is a schematic drawing of 
the ^-TPC. A /i-PIC with a detection area of 31 x 31 cm^ (TOSHIBA/DNP, 
SN060222-3) was used as the main gas-multiplier and as readout electrodes. A 
gas electron multiplier (GEM [27]) with an effective area of 23 x 28 cm^, which 
was manufactured by a Japanese company (Scienergy Co. Ltd.), was used as a 
pre-amphfier. The drift length was 31 cm and the detection volume was 23 x 28 
X 31 cm^. The electric field of the drift volume was formed with a field-shaping 
pattern on a fluorocarbon circuit board surrounding the detection volume. The 
detector was placed in a 3 mm-thick stainless-steel vacuum chamber filled with 
150 Torr CF4. The data-acquisition system, which is described in Ref. [28], was 
triggered by the signals of the /x-TPC itself. Hence, the absolute Z position was 
not measured, but only the relative positions. Typical operation parameters 
were V^_pic = 490V, Vgemb = 630V, Vgemt = 915V, and VDRiFT=8.25kV. 
V^-pic is the voltage supplied to the /i-PlC, Vgemt is the voltage supplied 
to the top of the GEM, Vgemb is the voltage supplied to the bottom of the 
GEM, and Vdrift is the voltage supplied to the drift plane. These parameters 
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were optimized to realize stable operation at a combined (/i-PIC x GEM) gas 
gain of 2300, which is high enough to maintain the detection efficiency of the 
recoil nucleus. We still suffer from discharge problems when operating a /i-PlG 
alone. 



2.2 Energy calibration and energy resolution 

We calibrated and monitored the energy with alpha particles generated by 
the ^°B(n, a)'''Li (Q=2.7MeV) reaction [23]. We set a glass plate coated with 
a thin 0.6 fim ^^B layer in the fi-TPC, and irradiated the /i-TPC with ther- 
malized neutrons. Alpha particles (5.6 MeV, 6.1 MeV, and 7.2 MeV) from the 
decays of the radon progeny were also used. The track-length and the energy- 
deposition correlation down to the DM energy range agreed with a calculation 
by SRIM[29]. 

The measured energy resolution in 5-8MeV was 50%(FWHM), which was 
dominated by the gain inhomogeneity of the /i-PIC. This would be the worst 
value when the gain inhomogeneity of the /x-PIC dominates the energy reso- 
lution. The energy resolution could be dominated by the the statistics of the 
primary electrons. In this case, the energy resolution at 100 keV with a CF4 
gas was estimated from a measurement with an Ar-C2H6 gas mixture, because 
it is not very easily measured due to a small photo- absorption cross-section 
and a long range of the electron. The energy resolution at 60 keV, measured 
with an Ar-C2H6 gas mixture, was 60% (FWHM)[26]. We caliculataed with 
AEcF, = A^Ar • i ito/wc^ y^^- ^^CF4 and A^Ar are the energy reso- 

lution at 100 keV with a CF4 gas and that at 60 keV with an Ar-C2He gas 
mixture, respectively. WcF4=54eV and WAr=26eV are the energies needed to 
created an ion-electron pair in a CF4 gas and an Ar-C2H6 gas mixture, re- 
spectively. Thus, the estimated energy resolution was 70% (FWHM) at 100 
keV with a CF4 gas, if the statistics of the primary electron dominates the 
energy resolution. The energy resolution in the DM energy range was therefore 
estimated to be better than 70%. 



2.3 Nuclear recoil detection efficiency 

The nuclear-recoil detection efficiency was measured by irradiating the /x-TPC 
with fast neutrons from a ^^^Cf source of 0.72 MBq placed at several posti- 
tions. We selected the tracks within the fiducial volume (21.5 x 22 x 31 
cm^). We also required at least three digital hit points for a track. A typlical 
track-length and energy-deposition correlation is shown in the left panel of 
Figure 2. A nuclear recoil band is seen along the X-axis. We set the upper 
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limit of the "nuclear recoil band" at 1 cm. We compared the spectrum after 
the selection with a simulated one, and the calculated ratio was the overall 
(detection and the selection) efficiency. The efficiency curve is shown in the 
right panel of Figure 2. The efficiency was about 40% at 100 keV, and the 
efficiency dependence on the recoil energy was phenomenologicaly fitted with 
an error function, e. The best-fit function was e = 1.0 • erf{{E — 45.8)/165.2), 
where erf{x) — Jq exp{—x^)dx is the error function and E is the energy 
in keV. 



2.4 Gamma-ray rejection factor 

We demonstrated that the track-length and the energy-deposition correla- 
tion measured with a gaseous TPC provided a strong gamma-ray rejection in 
our previous work[24]. We performed quantitive measurements with a large- 
volume /i-TPC for a WIMP-search measurement in a surface laboratory. We 
irradiated the //.-TPC with gamma-rays from a ^^^Cs source of 0.80 MBq, 
placed at 25 cm from the center of the detection volume. The measured track- 
length and energy-deposition correlations with and without the gamma-ray 
source are shown in the left and the right panels of Figure 3, respectively. 
The live time was 0.44 days for both measurements. A band along the Y-axis 
labeled "electron recoil band" is seen only in the results with a gamma-ray 
source, which is obviously due to tracks of the recoil electrons. Another band 
along the X-axis, seen in both results, is due to the tracks of the carbon and 
fluorine nuclei by the neutron background, and are labbeled "nuclear recoil 
band" . The spectra of these two measurements are shown in the left panel of 
Figure 4. The spectrum of the no-source run was subtracted from the that 
of the ^^''Cs run. The obtained "subtracted" spectrum is shown in the right 
panel of Fig 4. The subtracted spectrum is consistent with zero in the DM 
energy range, within the statistical error. We then derived the gamma-ray 
rejection factor by dividing the subtracted spectrum by a simulated electron 
recoil specrum. The measured gamma-ray rejection factor in the DM energy 
ranege was (0.57 ± 0.69) x 10^^. The error is statistical, which is dominated 
by the background in the surface laboratory. The 90% C.L. upper limit of the 
gamma-ray rejection factor was 1.5 x 10~^. 



2.5 Direction- dependent detection efficiency 

The yU-PIC has an anisotropic response with regard to the three-dimensional 
track directions. This non-uniform response is due to a difference in the de- 
tection method of the Z-axis (time information) and the other two axes (2- 
dimensional image taken by the /x-PIC). We measured this detection efficiency 
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dependence on the track direction, or a "response map", by irradiating the /x- 
TPC with fast neutrons from several positions that generated uniform recoils. 
We fitted the digital hit points with a straight line, and determined the direc- 
tions of each nuclear recoil track. Because we did not measure the track senses 
in this measurement, we fixed the sign of the x component to positive and 
determined the (x,y,z) vector of the track in the coordinates shown in Fig 1. 
The measured response map is shown in Fig 5. The map was normalized by a 
mean value. It can be seen that the correction factor is quite large because the 
fabrication technology of the large-sized //-PIC is still in an R&D phase, and 
the gain inhomogeniety is still large. The detector showed a relatively high 
efficiency for tracks in the X-Z and Y-Z planes, but showed a relatively low 
efficiency for tracks in the X-Y plane. 

We have not overcome intrinsic weak points of the TPC that the efficiency to 
the tracks parallel to the readout plane is relatively low. It should be noted 

that the absolute detection efficiency was measured by a method described in 
the pervious sections, and this response map was used to correct the relative 
efficiency. 

2.6 Position resolution and angular resolution 

We measured the position resolution using the tracks of the alpha particles 
measured in the energy calibration. We evaluated the position resolution using 
the same method as described in our previous work [26]. The measured posi- 
tion resolution was 0.8 mm (rms). The angular resolution of the carbon and 
fluorine nuclei in the DM energy range was estimated in the following manner. 
We calculated the samphng pitch of the carbon and fluorine tracks from real 
data. We then generated the hit points by a simulation based on the mea- 
sured position resolution and the sampling pitch. Simulated hit points were 
fitted with a straight line, and the angle between the fitted line and the real 
direction was calculated. The distribution of the calculated angle was fitted 
with a Lorentzian and 7 = 25°, where 7 is a parameter of a Lorentzian that 
corresponds to HWHM. We used this value for a direction-sensitive analysis. 



3 Measurement 

A dark matter search measurement was performed in a surface laboratory 
at Kyoto University {Lat. 35°2' N, Long. 135°47' E, third floor in a five- 
story building). The measurement was performed from November 1st, 2006 to 
November 27th, 2006. We set the /i-PIC plane horizontally and re-aligned the 
X- and Y-axis of the /x-TPC twice (three different indicated in Figure 1) 



6 
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-5° 


4.37 


40° 


5.25 


18° 


4.95 


18° 


i^^Cs 2.14 



Table 1 

Properties of the dark matter measurements. 

so as to cancel any remaining static errors. Table 1 lists the measurement 
properties. The total live time was 16.71 days and the exposure was 0.151 kg- 
days. The sum of the calibration and maintenance time was 7.51 days, and 
the dead time due to the data acquisition was 1.0 days. 

The yU-TPC was filled with 150 Torr CF4, and the vacuum chamber was sealed 
without any gas circulation.. The drift velocity, total gas gain, and the energy 
resolution were monitored every two or three days during the measurement. 
The drift velocity was 7.6 cm//xs at the beginning of the measurement and 
decreased to 6.9 cm/fis at the end. This decrease is thought to be due to the 
out-gas from the detector components. We used the measured velocities for 
the analysis. The mean gas gain in the entire detection area was 2300, and a 
time-dependent variation was not observed within a measurement error of 4% 
(rms). The energy resolution was also stable within a measurement error of 
7% (rms). The time-dependence was small compared to the energy resolution, 
itself, and the angular resolution in this measurement. 



4 Results 

The spectra did not show a statistically significant difference among the four 
conditions due to the strong gamma-ray rejection. We therefore combined the 
spectra obtained in the four measurements, including even the spectrum ob- 
tained in the gamma-ray rejection measurement with a ^^^Cs source. Figure 6 
shows the obtained spectrum in the DM energy range. The spectrum can be 
explained by a typical neutron ffux in a surface laboratory of 0(1 x lO^^n-cm^^- 
s~^). We initially derived the limits of the WIMP-proton spin- dependent cross 
section only from the spectrum (conventional method) . We used the parame- 
ters given in Table 2, and followed Ref.[30] for the calculation. The thin-dotted 
line, which is indicated by "spectrum only", in Figure 9 shows the exclusion 
limits. 

We next performed a direction-sensitive analysis. We created a sky map of 
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A2J(J+1) = 0.647 



Tabic 2 

Astrophysical and nuclear parameters used to derive the WIMP-proton limits. 

the north sky with the direction of the three-dimensional nuclear recoil tracks. 
Because we did not detect the sense of each track, the map was restricted to 
half of the sky larger and the southern part of the sky was folded. We plotted 
every nuclear event (1686 events in a 0.151 kg-days exposure) in Figure 7. 
The larger mark, labled "DM expected", is the direction of the solar motion 
or the Cygnus from which the WIMP-wind is expected. The WIMP-wind 
direction has a diurnal motion around the polestar. We calculated 6, the angle 
between the WIMP-wind direction and the recoil direction for each event, 
and drew the the cos 6' distribution shown in Figure 8. We compared this 
measured cos^ distribution with the expected WIMP-signal. The expected 
WIMP-signal was prepared by the following method. We first created a cos 9 
distribution for a given WIMP mass and a recoil energy while considering 
the energy resolution. The DM energy range was divided into 15 bins, and 
the energy and angular resolution of the detector were taken into account. 
We used 70% (FWHM) for the energy resolution and an angular resolution of 
7 = 25°. We followed Ref. [16] for the angular distribution and Ref.[30] for the 
count rate. We then created a sky-map expected by the WIMP events with 
this cos^^ distribution while considering the derection-dependence(Figurc 5) 
and the angular resolution. We used this sky map to create the expected 
cos 9 distribution for a given WIMP mass and energy bin. Figure 8 shows one 
of these cos^ distributions. Here, the data was fitted with the distribution 
for Af^=100GeV and 100-120 kcV bin, and the cross section that gave the 
minimum was 1.36 x lO'^pb and x^/d.o.f = 22.4/9. This best-fitted WIMP 
signal was rejected at the 90% confidence level by tests. The cross section 
that would give the minimum was calculated for each energy bin, and 
the smallest cross section was taken as the hmit for the given WIMP mass. 
The upper and lower panels in Figure 9 show the cross section limits and the 
corresponding values as a function of the WIMP mass, respectively. The 
limits are shown with a thick-solid line (labeled " direction-sensitive" ) and the 
values are shown by fiUed-square marks. The values at the 90% C.L. 
upper limit are shown by a dotted line in the lower panel in Figure 9. The 
WIMP signals were rejected at the 90% confidence level because the best-fit 



8 



values were all above the 90%-%^ line. The data was also fitted with a 
fiat cos^^ distribution and x^/d.o.f = 9.9/9 was obtained independent of the 
WIMP mass. The flat 0086* was not rejected at the 90% confldence level. It 
was demonstrated that the direction-sensitive method gave a limit slightly 
better than the conventional spectrum-only method. The difference between 
this result and other experiments is explained by the neutron background 
in the surface laboratory. This is the first limit set by a dark matter search 
experiment with a gaseous tracking device. As for the spin-dependent limit, 
this is the first result by a direction-sensitive experiment. It should be stressed 
that the dark matter search experiment with a gaseous tracking device can 
actually set a limit on the WIMPs. These results are very convincing and 
encouraging towards developing a future WIMP-wind detection with large- 
volume detectors. 



5 Discussions 



We have shown that WIMP-proton cross section limits can be obtained by 
a direction-sensitive method. However, much work is necessary to reach the 
sensitivity achieved by other experiments performed with solid and liquid de- 
tectors. The sensistivity is restricted by the neutron background in a surface 
run, therefore, we are going to perform a dark matter search experiment in 
an underground laboratory. The fast neutron fiux in Kamioka Observatory 
((1.9 ± 0.21) xlO^^n • cm~^ • s~^) is more than 3 orders of magnitude smaller 
than a typical fast neutron fiux in a surface laboratory. It it expected that 
the background sources, such as radioactive isotopes in the detector compo- 
nents, would limit the sensitivity in the first stage of an underground run. 
We are going to replace the detector components with low-activity materi- 
als. The angular resolution is not good enough to take full advantage of the 
direction-sensitive method over the "spestrum only" method. We are going to 
investigate the detector response using a neutron source of a known energy. 
We plan to study the operation of the //-TPC at a lower pressure of 30 Torr to 
lower the energy threshold. We are going to develop a larger readout area by 
placing many pieces of the of /x-PICs of 31 x 31 cm^. After these technological 
breakthroughs are achieved, we expect to reach the SUSY predictions with an 
exposure of more than one hundred kg-days, which we have calculated in our 
previous work[22]. 



9 



6 Conclusion 



We developed a three-dimensional tracking device, and performed a direction- 
sensitive dark matter search measurement in a surface laboratory. By using 
150 Torr carbon-tetrafiuoride (CF4) gas, we set the first limit on the spin- 
dependent WlMP-proton cross section by a direction-sensitive method. Thus, 
we have demonstrated that a dark matter search experiment with a gaseous 
tracking device can actually set a limit on WIMPs. Furthermore, we have 
demonstrated that this method will potentially play a certain role in revealing 
the nature of dark matter when a low-background large-volume detector is 
developed. 
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Fig. 1. Schematic drawings of the /i-TPC. /i-PIC and the GEM are set horizontally. 




Fig. 2. Measured track- length and energy-deposition correlations with a neutron 
source of ^^^Cf (left) and the detection efficiency dependence on the recoil energy 
(right). 
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Fig. 3. Measured track- length and energy-deposition correlations with (left) and 
without (right) the gamma-ray source of ^^^Cs. Two bands along the Y- and X-axises 
are seen in the ^^"^Cs run, while the electron recoil band is not seen in the no-source 
run. Several events out of both bands are thought to be the tracks of protons from 
the wall of the drift volume. 
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Fig. 4. Spectra taken with (thick) and without (hatched) the gamma-ray source 
^^''Cs (left) and the subtracted spectrum (right). 
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Fig. 5. Measured direction-dependent efficiency of the /x-TPC. X means the direction 
along the X-axis. The map is normahzed by the mean efficiency. 
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Fig. 7. Obtained sky map of the north sky. The recoil directions of the fluorine and 
carbon nuclei are used. A raw map and a corrected map by the relative efficiency 
dependent on the directions are shown. The thick line indicates the direction of the 
solar motion (direction of the WIMP-wind). 



17 



400 



350 



300 



250 



200 



150 



100 



50 



MylOOGeV (100-120keV) 



13600pb 
X =22.4 



simulation 
data 



0.2 



0.4 0.6 
cos 



0.8 



Fig. 8. Measured (with error bars) and expected (histogram) distribution of the 
angle between the recoil direction and the WIMP direction. The expected histogram 
is that with M^=100GeV, 100-120 keV bin, and 1.36xl0^pb. 



18 



SD 90% C.L. upper limits 



10° 
10^ 
^ 10^ 

- 10"* 



Q. 
Q- 



1 lo-^ 

^ 10^ 
10^ 



1Q 



■1 
'30 
20 
10 




'\ This work (CF4 


spectrum only) ■ 


This work (CF4 


direction-sensitive) 


: DAMA(Xe) 






CDMS(Ge): 


^^Tokyo 2006(CaF2) ^^^^^j^^j^^^^^^^^^^a^t^-- 


\ DAMA(Nal)..--': , 




" ■ "Bm 

90% limit ■ 


WIMP 

■ ■ ■ ■ ■ ■ II 


• • ••••• • 

1 


• • • • • • 

flat BG 



10' 



M 



10-= 

WIMP 



[GeV] 



10' 



Fig. 9. 90% C.L. upper limits on the WIMP-proton spin-dependent cross section 
(upper) and values (lower) as functions of the WIMP mass. The thick solid 
and dotted lines show the limits obtained with and without the direction informa- 
tion, respectively. Limits from other experiments ( DAMA(Xe)[3], DAMA(NaI)[4], 
NAIAD[5], Tokyo CaF2[7], SIMPLE 2005[10], PICASS0[11], and CDMS[12]) are 
shown for comparison. The filled squares show the minimum values of the best- 
fit WIMP distribution, the filled-circles show the best-fit flat cos 9 distribution, and 
the dotted line show the values at the 90% C.L. upper limit. 
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